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Abstract 

The goal of this paper is to give a detailed analytical description of the global dynamics of 
A'' points interacting through the singular logarithmic potential and subject to the following 
, symmetry constraint: at each instant they form an orbit of the dihedral group Di of order 

■ 21. The main device in order to achieve our results is a technique very popular in Celestial 

^ ' Mechanics, usually referred to as McGehee transformation. After performing this change of 

, coordinates that regularizes the total collision, we study the rest-points of the flow, the invariant 

' manifolds and we derive interesting information about the global dynamics for I — 2. We 

I observe that our problem is equivalent to studying the geometry of stationary configurations 

00 . of nearly-parallel vortex filaments in three dimensions in the LIA approximation. 
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■ Introduction 



Equations of motion for interacting point vortices were introduced by Helmholtz in a seminal 
paper published in 1858. Towards the end of the paper he introduced the point vortex model, 
by considering the trace of the point of intersection of a family of infinitely thin, straight parallel 
Q>^ ' vortex filaments with a plane perpendicular to one (and then to all) of them. A large and still 

OO . growing body of literature has focused on the study of point vortices, seen as vorticity monopoles 

' in a two-dimensional ideal fluid. One may think of them as playing a role in ideal hydrodynamics 

similar to that played by point masses in Celestial Mechanics. For a review see the book by Newton 
' [NewOl] and references therein. 

A comparatively much smaller body of literature has been devoted to the study of vortex 
filaments in a three-dimensional ideal fluid, a subject probably closer to Helmholtz's original ideas. 
Based upon the equations governing the evolution of vorticity in three dimensions. Da Rios derived 
in 1906 the localized self-induction approximation (LIA) describing the approximate motion of an 
, isolated filament and later re-derived by Arms and Hama in 1965. In 1972 Hasimoto introduced 

^ ' a change of coordinates which takes the familiar Frenet-Serret formulas for the geometry of a 

5^ , curved filament and the localized self-induction equation governing its dynamics into the nonlinear 

Schrodinger equation, a completely integrable infinite dimensional Hamiltonian system. Building 
on these results, Klein, Majda and Damodaran in [KlMaDa95] derived in 1995 a simplified model 
describing the time evolution of A^-vortex filaments nearly but not perfectly parallel to the z- 
axis. According to this model the motion of A'^-interacting nearly parallel filaments is given by the 
following coupled partial differential equations: 



(1) -dt^, =^jdl^j+J2 

* Dipartimento di Matematica, Universita del Salento, Italy, email: francesco.paparellaaunisalento.it. 

t Dipartimento di Matematica, Universita del Salento, Italy, email: alessandro .portaluriOunisalento . it. 
Work partially supported by"Progetto 5 per mille per la ricerca" (Bando 2011). "CoUisioni fra vortici puntiformi e 
fra filamenti di vorticita': singolarita', trasporto e caos.". Work partially supported by the PRIN2009 grant "Critical 
Point Theory and Pcrturbative Methods for Nonlinear Differential Equations" . 



1 



for N := {1,...,A^}, Tj e M* and where each '^j is a complex- valued function. Among all the 
solutions of the system (1) a special role is played by the stationary solutions. 

In this paper we are interested in studying the geometry of the stationary solutions of (1). 
Therefore, by substituting in (1) 'i>j{a,t) with qj{a), where q : R — > we reduce (1) to the 
following system of coupled ordinary differential equations for stationary vortices filaments: 



where ' denotes differentiation with respect to the arc-length parameter cr. By multiplying each 
equation in (2) by Tj and defining the potential function U as 



where q :— (qi, . . . , qn), the ODE in (2) can be written as 



Here F is the diagonal block matrix defined by F = [Tij] and Tij = r^dijl2 where I2 denotes 
the two by two identity matrix. By interpreting the parameter cr as a time-like coordinate, the 
equation (4) can be seen as describing the dynamics of point masses on the plane, interacting with 
a logarithmic central potential. Therefore, from now on, we shall refer to the parameter cr as a 

time-parameter. Let F* — X^j^i^j- We define the center of vorticity C as C := Y^'l2jLi^'j^j- 
From (4) it is straightforward to check that the center of vorticity satisfies the conservation law 
C = 0. A severe difficulty in investigating a system such as (4) is due to the presence of singularities 
both partial and total: physically they represent collisions between some or all of the vortices. This 
problem motivates the introduction of a change of coordinates that regularizes the equation of 
motion. A change of coordinates of this sort was recently used by Stoica & Font in [StFoOS] for a 
singular central log-potential problem which arises in galactic astrodynamics. This is an example 
of McGehee transformations ^ a regularizing change of variables currently popular in the field of 
Celestial Mechanics and first introduced in 1974 by McGehee [McG74] in order to study orbits 
passing close to the total collapse in the coUinear three-body problem in R^. 

McGehee transformations consists of a polar-type change of coordinates in the configuration 
space, a suitable rescaling of the momentum and a time scaling. The idea behind this non Hamil- 
tonian change of coordinates is to blown up the total collision to an invariant manifold called 
total collision manifold over which the flow extends smoothly. Furthermore, each hypersurface of 
constant energy has this manifold as a boundary. The effect of rescaling time, is to study some 
qualitative properties of the solutions close to total collision. In fact, by looking at the transforma- 
tion defined in equation (9), it readily follows that the effect of this transformation is to slow the 
motion in the neighborhood of the total collapse, which is reached in the new time asymptotically. 

However, we point out that there are several crucial differences between the classical McGehee 
transformations in use in celestial mechanics and those appropriate for our case. The most important 
one is related to the lack of homogeneity of the logarithmic nonlinearity. This breaks down some 
nice and useful properties of the transformation. For instance it is not possible as in the iV— body 
problem, to recover the global dynamics of orbits passing arbitrarily close to total collapse by merely 
looking at a suitable Lyapunov function defined on the total collision manifolds whose only rest 
points represent the central configurations of the bodies in gravitational interaction. (Compare 
[McG74], [Dev80] [Dev81], for further details). 

Nevertheless, it is still possible to regularize the vector field and therefore we can still carry out a 
detailed analytical description of the rest points, of the invariant manifolds and of the heteroclinics. 
The ability of these coordinates to give insight on the dynamical properties of the problem becomes 
apparent when studying a class of symmetric solutions that we call dihedral equivariant (Section 2). 
In the new coordinates we were able to investigate the global dynamics of the problem and to prove 
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some important features which were not at ah obvious when stating the problem in conventional 
Cartesian coordinates. For example, we could prove that a dihedral equivariant configuration of 
four vortices is always bounded (Theorem 5.2), even if the standard argument based on total energy 
conservation is unable to rule out unbounded solutions, as it does for central potentials like those 
studied in [StFoOS]. 

Another important feature is the fact that every dihedral equivariant solution experiences a 
collision within a finite time (Theorem 5.3). This may be either a total collapse (simultaneous 
collision of all the vortices) or a binary collision (simultaneous pairwise collisions of the vortices). 
Only the second type of solutions is generic, in the sense that it is generated by a set of initial 
conditions having full Lcbesgue measure (Theorem 6.6). By using a recent result proven in [CaTe] 
we were able to study the dynamics for arbitrarily long times, by defining generalized solutions that 
continue by transmission after a binary collision. 

On the other hand, an interesting qualitative feature that the present problem has in common 
with the gravitational n-body problem, is that the total collapse can only be reached in central 
configuration. In Celestial Mechanics this is a well-known fact first proven by Sundman for the 
three body problem and generalized by Wintner in the general case. Recently, the authors of 
[BaFeTeOS] proved that this result also holds for a class of weak potentials, including the logarithm. 
We recover this property as a dynamical feature of the equations written in McGehee coordinates. 
It is worth noticing that for genuine point vortices (and not for stationary configurations of vortex 
filaments, as in our case) total collapse need not happen along central configurations (that is, with 
self-similar motion; see [NewOl] for further details). 

Finally, our results hint at the presence of a chaotic dynamics of the generalized solution. A 
challenging problem would be to investigate the symbolic dynamics or to say something about the 
topological or geodesic entropy after a global regularization of the singularities. The dynamics of 
the non stationary solutions of the PDEs (1) also remains to be studied, even for initial conditions 
close to the stationary solutions investigated in the present work. 
The paper is organized as follows: 
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1 McGehee coordinates and regularization 

In this section we develop the McGehee-like transformations that we need in order to study the 
dynamics close to total collapse. In the first subsection, we describe the general framework and we 
fix notations. Then, in the following subsection we repeatedly perform coordinate changes, and we 
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re-parameterize time, until we arrive at a set of equation of motions which is suitable to study the 
dynamics of our problem. 



1.1 General set-up and McGehee coordinates 

Let > 2 be an integer. Let denote the origin in and let Fi, . . . , Fjv be N positive numbers 
(which can be thought as strength of the vortex filaments) . The conservation law of the center of 
vorticity implies that there is always an inertial reference frame where the position of the center of 
vorticity is at the origin. In this reference frame we can identify the configuration space Q for the 
problem with the subspace of M?^ defined by 

N 

Q:={g=(gi,...,qjv)eM2A^: ^F.g, = 0}. 

For each pair of indexes i, j € N let Aj^ denote the collision set of the i-th and j-th vortex filament; 
namely Aij = {</ G : Qi = qj}. We call collision set the subset of the configuration space given 
by A := [Ji^j A^j and reduced configuration space, the set Q Q \ A. The Newton's equations 
in Hamiltonian form can be written as: 



(5) 



rq = p 

. _ mj_ 

^ dq 



where the Hamiltonian function H is defined by: 

(6) ff:r*(Q)^R H{q,p):=]^{T-^p,p)-U{q), 

on the phase space T*{Q) (the cotangent bundle over the configuration space). Explicitly T*{Q) = 
Qx P for 

N 

P:={p=(pi,..-,Piv)GR''^| Ep^^O}- 

The differential equations in (5) then determine a vector field with singularities on R^^ x R^^, or a 
real analytic vector field without singularities on (R^^ \ ^) ^ R^^. The vector field given by (5) is 
everywhere tangent to Qx P and so this 4(A'^— 1) dimensional linear subspace is invariant under the 
fiow. We henceforth restrict our attention to the fiow on the phase space Qx P. Consequently H is 
an integral of the system. This means that the level sets Ylh '■= H~^{h) r){Q x P) are also invariant 
under the flow (5). We observe that E/j is a real analytic submanifold of Q x P having dimension 
4(A^ — 1) — 1. The flow, however, is not complete. In fact certain solutions run off in finite time. 
This happens in correspondence to any initial condition leading to a collision between two or more 
vortex filaments: the corresponding solution meet A in finite time. We shall call total collapse, or 
total collision the simultaneous collision of all the vortices. Because the center of vorticity has been 
fixed at the origin, if a total collapse happens, it must occur at the origin of Q. 

We introduce the angular momentum. Denoting by (R^,a;o) the two dimensional symplectic 
vector space, the angular momentum is given in terms of the standard symplectic structure as 

N 

(7) M ■.= Y,T-'u:o{q„P,). 

We observe that if qj := (<7j,0) £ R'', pj := {Pj,0) <= R'' then the angular momentum defined in 

(7) agrees with 

N 

(8) Meg = ^FTi(q^. xp,). 
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An important class of solutions for many problems of celestial mechanics is called s elf- similar . 
This class is characterized by the fact that up to rotation and homothecy the configuration of all 
point masses is constant in time. Self-similar solutions exist also for our problem. Before proceeding 
further, let us first define a central configuration. 

Definition 1.1 A point s e M^^, such that {Ts,s) = 1 is called central configuration if 



for some fi £ 



We now characterize the class of self-similar solutions. 

Lemma 1.2 Let s be a central configuration, and p be a positive smooth function such that 



Then q{(T) = p{cr) s is a solution of (4). 
Proof. By a direct substitution, we have 

Tq = Tps = T^,p-'s = p-' — {s) = —{q\ 

where the last equality follows by the fact that the gradient of [/ is a homogeneous function of 
degree —1. □ 
Let fi,(p2 be two monotonically increasing, unbounded, smooth functions on the positive half- 
line, vanishing when r — 0+. We define the following McGehee-like coordinates: 



(9) 



WqW = fiir) 

s = ql \\q\\ 



The vector s lies on the 2{N — 1) — 1-dimensional unit round sphere centered at the origin, that we 
shall call the shape sphere S. As a direct consequence of the fact that the center of vorticity is fixed 
at the origin, it follows that a total collapse may only happen at r = 0. This constraint defines the 
collision manifold. 



1.2 Equation of motions in {r,s,z)- coordinates 

In order to write the equations of motion using the McGehee coordinates, it will be useful to define 
V := {T~^z,s) and to rescale time according to [(pi(p2\~^da = dr. In these new coordinates the 
submanifold of constant energy h is given by 

(10) E,. = {(r,s,2) gM; X 5x R2 I ^ 2^l{r)[h + UiMr)s)]} . 

By taking into account the expression (3) for the potential, we have 



<P2 



{r)U{Mr)s) := -^l{r) ^ log H = 



47r 
47r 



(logVlW +l0g|lSj - Sfcll) 



In order S/i to be not empty we choose the functions Lpj as follows: 

ipi{r) := 7y(r)e-['^^('"'l"', r/(r) £ 0{r^), k > 0. 
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With this choice, we get 1 := hni^^Q+ —ip2{r)\og(pi{r) and by this it follows that the stratum E;, 
meets the collision manifold r = along the submanifold 



A := <{ (0, s, 2) e R; X 5 X R2 I ^{T^^z, z) 



G 



where G := J2k 



r,r 



We observe that A does not depend on the energy level h. Therefore, all 



strata T,h share the same boundary at the collision manifold. Following the authors in [StFoOS] 
one possible choice is the following: 



(11) 




With this choice, the Hamiltonian system (4) becomes the following system of ordinary differential 
equations 



(12) 



r-'^z- vs 



1 9U , , 

:vz + T^(s) 



dq 



where, in the last equation, we have used the identity re dqU{q) ~ dqU{s) which follows from 
the homogeneity of degree —1 of dqU , together with (9) and (11). 



2 The dihedral problem 

Let us recall some basic facts about the dihedral group, seen as a map of R'^ into itself. For further 
details we shall refer to [FePoOS]. Let R'' = C x R be endowed with coordinates {z,y), z G C, y £ R. 
For / > 1, let Q denote the primitive root of unity Q — g^'''/'; the dihedral group Di C 5*0(3) is 
the group of order 21 generated by 

Ci- {z,y) {Ciz,y), k: {z,y) 1-^ iz,-y), 

where z is the complex conjugate of z. The non-trivial elements of Di = {Q , k) are the I — 1 rotations 
around the Z-gonal axis , j = 1, . . . , Z — 1 and the I rotations of angle tt around the / digonal axes 
orthogonal to the /-gonal axis (see figure 1) k, j ~ I, . . . ,1. In figures 1(a) and 1(b) we show the 
upper-halves of the fundamental domains for the action of Di restricted on the unit sphere for I = 2 
and 1 = 3, respectively. 

The action of Di, restricted on the fixed subspace (R"*') ' = R^, generates a dihedral configu- 
ration of vortices for any given Qq G M.'^. If we assign the same circulation to all vortices, then we 
can express the logarithmic potential acting on a single vortex as 

(13) uiqo) = - logdko-.ggoll), 

aeDi-^{i} 

where, without further loss of generality, we have taken Ff = 27r/L A solution of equations (5) such 
that the vortices lie on a dihedral configuration at all times, is a dihedral equivariant orbit. 

Lemma 2.1 The angular momentum of any dihedral equivariant orbit is zero. 

Proof. Using (8) and Hamilton's equations (5), it follows that the angular momentum is: 

N 
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(a) 1 = 2 (h) 1 = 3 

Figure 1: Dihedral groups Di, with the upper half of the fundamental domains in white. 

By hypothesis, at each instant the configuration of the vortices is given by the action of the dihedral 
group Di; Qj = g.jq for gj G Di. Hence 

qj X qtj g^q x g^q gj{q x q) . 

We observe that q x q = (0, 0, y), for some real number y. Therefore, using (2), we have 



M = Y, g, {qy.q)= /(0, 0, y) + Z(0, 0, -y) = 0. 

i=i 



This conclude the proof. 



□ 



Remark 2.2 As consequence of the Sundman-type estimates proved by authors in [BaFeTeOS], for 
a large class of potentials including the logarithmic one, a necessary condition in order to have the 
total collapse is that the total angular momentum of the system should be zero. Therefore as direct 
consequence of Lemma 2.1, we can conclude that total collision orbits can occur. 

2.1 The logarithmic dihedral potential 

For z — pe'", with p > and aj = 27rj//, j = 1, — 1, it is easy to verify that: 

\z — z\ = 2p |sina| 
|z-C/z| = 2p|sin(^) 

\z-Cfz\ = 2p smi^-^ 

Therefore, for a G (0, tt//) (that is, for q in the fundamental domain) the potential (13) becomes: 



(14) 



U{p,a) :=-log((2p)2'-i)-lo, 



i-i 



sin a Y[ sin ( y) sin - aj 



Since 11^=1 sin \J~f) ^ 2i-i ^^'^ ^■^ taking into account the multiple-angle formula 

linx) = 2"-' Yl sin f—+x) , 

fc=o ^ " 



sm 
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we have 

i-i 

sin a JJ^ sin 

These calculations lead to the following definition. 



JIT 

T 



sin(Za) 



Definition 2.3 We define the dihedral logarithmic potential as : 

(15) U{q) :=-logp2'-i_log(2/sin(;a)), 
and its angular part as: 

(16) C/(a) -log(2;sin(;a)). 

Lemma 2.4 (Planar 2l-gon) The dihedral problem admits exactly one (up to permutation of the 
vortices) central configuration, which is given by the vertices (e^^'^^-'^^'^'/^^'), 0) of a regular 2l-gon. 

Proof. From the definition (1.1) it follows that central configuration corresponds to a critical point 
of the potential restricted to the shape sphere. From (16) it follows that 

U'ia) ^ 

sm(ia) 

7r 

whose critical point in the fundamental domain of the shape sphere ^® = From U"{ac) ~ P 

it also follows that the point Uc is a minimum. □ 
We also observe that 

1. U{ac) < for Z = 2,3; 

2. U{ac) > for ? > 4. 
By using the identity 



sin(na) = ( j cos*"' a sin" ^ a sin i-{n— k)7r 



k=0 



the expression (16) may be written in terms of the local parameterization of the unit sphere s{a) 
(SI7S2) := (cos a, sin a) as follows 



(17) 



[/(s) = -log 



s\s2 sm 



in(f(^-.)) 



For j G 2, qj = Sj||q|| and hence the dihedral potential given in definition 2.3 becomes: 
(18) C/(q) = -log 



2.2 The geometry of the energy hypersurfaces 

From (10) let us define E{h, r, s) := 2r^{h + U{re^^^'' s)). More explicitly E is given by: 

E{h, r, s) := E{h, r) + 2r'^ U{s), 

where E is given by 

(19) E{h,r) :=2[/ir2 + (2/-l)(l-r2logr)] . 
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1.0 1.5 2.0 2.5 3.0 



Figure 2: Zero set of the function (r,a) i— >■ E{0,r,a) and different value of the total number of 
point vortex. 

Remark 2.5 The boundary of the regions where the motion occurs is given by E(h,r,a) = 0; more 
precisely, since the kinetic term {T~^z,z) is a positive definite quadratic form, for any fixed energy 
level h the motion is possible where: 

Eih,r,a) := 2 [/i + (2Z - 1)(1 - log r)] - 2r2 log(2/ sin(;a)) > 0. 

The shape of the curve E(h,r,a) = is key in order to understand the dynamics of our problem. 
Therefore the hypersurface corresponding to the energy level h is given by 

(20) T,h = |(r,a,2) e R; x (0,7r/0 x | {p-'^z^z) = E{h,r,a)^ . 

We also observe that S/i meets the boundary r ~ along a submanifold given by 

A := {iO,a,z) eR*+x {0,n/l) x | {T-^z,z) = 2(2/- 1)} . 

The set A is diffeomorphic to an open cylinder and it represents the component over the fundamental 
domain of the total collision manifold. As we shall see in the following, the total collision manifold 
is homeomorphic to a two dimensional torus. 




Figure 3: Energy hypersurfaces for different values of the energy level h and the zero plane. 

Let us introduce a further change of coordinates. In the dihedral problem the shape sphere re- 
duces to 5^, therefore it is natural to use the parameterization a I—)- s(a), where s(a) = (cosa,sina). 
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Following [StFoOS], we exploit the constraint (20) to parameterize the we also introduce a new local 
parameterization the momentum coordinates z with an angle in the following way 



(21) 



~ /27r , /27r . 

E\h^ r,a) ( W — cos ip, y — sm ^p). 



Note that with the choice of the strength of the circulations made in Section 2, we have F ^ = 1/2ttI. 
We also use a new time-like variable C, defined by 



dr = VEdC. 

With the above parameterizations, the system (12) reduces to 



(22) 



dr 
d( 

da 

dc 



27r r2 + 2 
I 



E{h,r,a) (cos a cos ip + sin a sin ip) 



E(h, r, a) (sin cos a — sin a cos "0) 



r^{dq^U{s{a)) cosip — dq^U{s{a)) sin V')- 



are 



The following result is an obvious consequence of equations (22). 

Lemma 2.6 The collision manifold at r = and the zero velocity manifold E{h,r,a) 
invariant manifolds. Moreover: 

1. On the collision manifold r — the dynamics is given by 

da I , , s . , s dib 

^=-V^2(2^-l)sm(0-a), -| = 0; 

2. On the zero-velocity manifold E(h, r, a) = 0, the dynamics is obtained by integrating the last 
equation in (22), keeping r, a as constants. 



3 Flow and invariant manifolds: the Klein group D2 

For I = 2 (the Klein group) we may carry out a detailed study of the global dynamic of the problem. 
We observe that the dihedral potential reduces to: 



whence it follows dqU{s) 



[/(g) = -log(8|lq|l(?i<j2) 
h si, h S2 I . In this case the equations of motion become: 

Sl 32 J 



(23) 



1 



dr 

dC " 7^r2 + 2 
1 



E{h, r, a) cos(?/' — a) 



da 

dC 

dip 
dC 



E{h, r, a) sm{i/j — a) 



-r 

7^ 



2 cos(-0 + a) 



sin('0 — a) 



sin 2a 

The rest points of (23) correspond to the solutions of the following systems: 



r = 

sm{tp — a) = 



r^O 
E^O 



E = 

2 cos{t/j + a) 
sin 2a 



sm{Tjj — a) 



However, it is readily seen that the second system has no solutions, as the conditions E = and 
r = are incompatible. With straightforward calculations we obtain the following result. 
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Figure 4: Rest points and invariant manifolds for the D2 problem. The green cylinder is the total 
collision manifold r = 0; the blue and the red lines lying on it are and 3^2: respectively. The 
four light grey manifolds are the surfaces defined by E{0,r,a) = (there is one for each of the 
four domains of the group D2); the red and the blue lines lying on them are =^3 and ^4 in the 
fundamental domain, and their images in the other domains. The blue/red color denotes families 
of rest points which are transversally linearly stable/unstable on their own invariant surface. See 
the text for the stability in directions transverse to the invariant surfaces. The four yellow planes 
correspond to binary collisions. 

Lemma 3.1 The equilibria of the vector field given in (23) consists of four curves, two belonging to 
the collision manifold and the other two on the zero- velocity manifold. In the coordinates {r,a,ij}), 
the curves on the collision manifold are given by 

(i) ^1 EE (0,a,a); 

(a) ^2 = (0,a, a + 7r) 

The curves on the zero-velocity manifold are given by 

(Hi) ^3 = (r,a, arctanm(a)); 

(iv) ^4= (r, a, arctanm(Q!) + tt), where 

, . cosa(sin^ a + 1) 

m[a) := -. 

sniQ;(cos'^ a + 1) 

and the pair (r, a) satisfies the equation E{h,r,a) = 0, which explicitly reads 

log(4sin2a) — h-\- 3(r^^ — logr). 

The spectrum and the eigenvectors of the Jacobian matrix of (23), evaluated at a rest point, 
gives useful information on the local dynamics close to the curves ^1, ^2, ^3, ^4- The logic 
of the calculation is elementary, but the computations occasionally become too large to be easily 
manageable. With the help of a computer algebra system (Maxima 5.21.1) wc arrive at the following 
results. 

Lemma 3.2 The curves of equilibria and ^2 (on the collision manifold) given in Lemma 3.1 
are degenerate. More precisely: 

1. at each point of the curve the spectrum is given by s{^i) — {Ai, A2} where Ai — 6/\/t^ 
and A2 = 0. Furthermore the algebraic multiplicity of Ai is h(Xi) — 1 and the algebraic 
multiplicity 0/ A2 is h{X2) ~ 2. The eigenvectors are 

Vl = (0,1,0) V2 = (1,0,0) -ua = (0,1,1) 
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where Vi is associated to Ai andv2, are associated to A2. The vector is tangent to 

2. At each point of the curve 3^2 the spectrum is given by 3(^2) = {Ai, A2} where Ai = — 
and A2 = 0. The multiplicities and the eigenvectors are the same as for 

Lemma 3.3 The curves of equilibria and ^4 (on the zero velocity manifold) given in Lemma 
3.1 are degenerate. More precisely: 

1. at each point of the curve iPy, the spectrum is given by three distinct simple eigenvalues, 
namely 3(^3) = {//i, /^2, /^s} where 

Ml =0 

[3 cos{^ - a)eS/'''+2'' + 64 cos(2a) cos(V' + a)r6] 



r^ [ cos(-0 — a) 6^/''^+'' + 8 sin(i/) + a) r^] 



g3/r2+;i VTT 



and 



£;o(r, a)=0, — =8111(7/1-0;). 

sm la 



The eigenvector associated to /is is 1x3 = (0,0,1), which lies on the invariant manifold 
E(h, r, a) = 0. The explicit expressions of the eigenvectors Ui and U2 associated to fii and 
to /i2 are not short, and we omit them. However, Ui is tangent to the curve ^3 and 1x2 is 
transverse to E{h, r, a) = 0. 

2. At each point of the curve ^4 the spectrum is 0(^3) = {/.ii, — yU2, — ^^3}. The eigenvectors are 
the same as for . 

Remark 3.4 For any a € (0,7r/2) it is fj.2 < and ^3 > 0. In fact, from the explicit expression 
of if,' {a) given in (3.1) we have 003(1/) — a) > 0, sin{ij) + a) > 0, cos(2a) 003(1/ + a) > 0. The last 
inequality stems from 

7r/4 < ?/'(q;) < 7r/2 - a if < a < 7r/4; 
7r/4 > > 7r/2 - a if 7r/4 > a > 7r/2. 

Remark 3.5 The above results show that the stability properties of the equilibria do not depend 
on the value of the total energy h. In other words, changing the parameter h does not lead to 
bifurcations. 

At each rest point let 113 denote by , respectively the (linearly) stable and unstable 
manifolds, and by the center manifold. The following results hold. 

Proposition 3.6 Any rest point in S^j, and ^4 is a degenerate saddle. More precisely we have 

dimM^"(^3) = 1, dimVK''(^3) = 1, dimiy°(^3) = 1. 

^. dimM^"(^4) = 1, diml^"(^4) = 1, dimM^0(^4) = 1- 

Proof. The statement is an immediate consequence of Lemma 3.2. Note that ^3 and ^4 are the 
center manifold of each of their own points, and that they are neutral, in the sense that each of 
their points is an equilibrium. □ 

Proposition 3.7 For each a, the two equilibria 

{0,a,a) e ^1 

and 

{0,a,a + TT) e 3^2 

are nonlinear degenerate saddles. 
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Table 1: Dimensions of the invariant manifolds for the equilibria belonging to the curves ^i, ,^2, 



(24) 



sin(-(/; — a) 



1. dimWi^i) ==1, dimM^O(^i) = 2. 

2. dim W {,^2) =1, dimW^0(^2) = 2. 
Proof. The equations of motion for an initial condition on the collision manifold reduce to 

dr 

da 
dC 

It is then evident that the unstable manifold of the equilibrium (0, ao, ao) G ^1 is the line (0, a, ao), 
and analogously (0, a, ao + tt) is the stable manifold for any point (0, ao, ao + tt) g ^2- 

The linear stability analysis of Lemma 3.2 shows that each point on ^1 and ^2 has a two- 
dimensional center subspace, and hence a two-dimensional center manifold. In fact I^i and ^2 are 
neutral submanifolds of the center manifold, in the sense that each of its points is an equilibrium. 

For brevity, here we shall not perform a formal center manifold reduction. In order to determine 
the stability of the dynamics on the center manifold we simply observe that the dynamics projected 
along the direction singled out by the eigenvalue V2 = (1,0,0), which is transverse to the total 
collision manifold, is given by the first equation in (23). From that, we have dr/dC, > in a 
neighborhood of any point of ^1 and dr/dC, < in a neighborhood of any point of i3^2- 

Therefore the equilibria of both and ^2 are nonlinear saddles. □ 



4 Heteroclinic connections and homothetic orbits 

The existence of heteroclinic connections on and between the invariant manifolds helps us to develop 
a global understanding of the flow. 

Lemma 4.1 The flow on the total collision manifold is totally degenerate. That is: 
(i) W^{Pi) = W%P2); 
(it) W^{P2) = W'{Pi); 

where Pi £ and P2 G ^2 o-fe chosen in such a way that the second coordinate of the two points 
is the same. 

Proof. The proof of this result follows by a straightforward integration of the equations of motion 
(24) valid on the total collision manifold r = 0. Note that the singularity of the potential at 
a = j7r/4, j = 0, • • • , 3 vanishes on the total collision manifold, therefore on the heteroclinic 
connection between Pi and P2, we have a G [0, 27r). □ 
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Figure 5: Total collision manifold and equilibrium curves. 



Remark 4.2 By taking into account the identification of the opposite edges of the rectangle {a,Tp) 
with the same orientation, the curves of equilibria can be identified with two closed curves on the 
torus and the heteroclinics are represented by arcs joining one point on a unstable (red) closed curve 
with the corresponding point on the stable (blue) curve as shown in Figure (5). 

Next we describe the flow on the zero velocity manifold E(h, r, a) = 0, where the equations of 
motion reduce to 



2cos('i/; + q) . , , , 
— sm^ip — a) 



sin 2a 

Lemma 4.3 The flow on the total collision manifold is totally degenerate. More precisely 

(i) 1^"(-P3) = ^^^(-Pi); 

(ii) W^{Pi) = W'{P:i); 

where P3 € ^3 and P4 € ^4 are chosen in such a way that the first two coordinates agree. 

Proof. From Lemma 3.1 we have P4 = P3 + (0,0, tt). From the equations of motion (25) follows 
that the flow lines on E{h,r,a) = are straight lines parallel to the i/j-axis. Therefore they are 
heteroclinic connections between each point P3 S ^3 and the corresponding P4 € ^4. □ 
By choosing (a,V') = (7r/4, tt + 7r/4) or {a,tp) = (7r/4, 7r/4) the last two equations in (23) are 
identically zero. This means that the constant pair (a,"*/") = (7r/4, tt + 7r/4) is a constant solution of 
the subsystem obtained by projecting the system of ode's onto the (a, '0)-plane. By summing up, 
the following result holds. 

Proposition 4.4 There exist two connecting orbits between the total collision and zero velocity 
manifold. More precisely 

1. The curve f7i3(C) = (^(C)j '''/4, 7r/4) is an heteroclinic joining and H^-}, for r solution of 

dr 1 r^ ~,, 

2. the curve 7724(C) = (''(O? ^/4, tt + 7r/4) is an heteroclinic joining ^4 and ^2 for r solution of 

dr 1 r^ 

where E{h,r,TT/A) = 2\hr'^ + - logr)] - 2r'^ log 4. 

Remark 4.5 From Lemma 2.^ with I ~ 2 it follows that the heteroclinic connections 7713 and 7724 
are self-similar orbits where the vortices are located at the vertices of a square. Their projection on 
the shape sphere is a central configuration. Moreover, these are the only self-similar orbits: from 
(23) it follows that no other solution curve has a(C) = 7''/4 when C ranges on an interval of non-zero 
length. 
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Figure 6: Heteroclinic connections among the points (0, 7r/4, 7r/4), (0, 7r/4, tt + 7r/4), (r, 7r/4, 7r/4), 
(f, 7r/4, TT + 7r/4) where f is such that E{h, f, 7r/4) = 0. 

Remark 4.6 A classical result of Sundman [Sun09] proved for the three-body problem in three 
dimensions that an orbit ending in triple collision asymptotically approaches a central configuration. 
The validity of Sundman-type asymptotic estimates for collision solutions is established in a recent 
paper [BaFeTeOS, Theorem 5, Example2] for a wide class of dynamical systems with singular forces, 
including the classical N-body problem, quasi-homogeneous and logarithmic potentials. Applying 
these results to our problem, we deduce that the only point of the total collision manifold which 
is (asymptotically) reachable by an orbit originating from an initial condition having r > is 
(r,a,ip) = (0, 7r/4, 57r/4), while (r^a^tp) = (0, 7r/4, 7r/4), is the only point reachable extending the 
orbit backward in time. In this sense, the collision manifold, with the exception of these two points, 
is dynamically disconnected from the region r > 0. 

Figure 6 shows graphycaUy the heterochnic connections between the invariant manifolds r = and 
E = 0. 

Remark 4.7 We stress that, while an initial condition on the heteroclinic rji^ is unable to complete 
the loop up to the total collision manifold (it just asymptotically reaches a fixed point on ^3), the 
dynamics expressed in Cartesian coordinates by (5) does perform the loop in a finite time. In 
fact, the points on the curves and ^4 of the zero velocity manifolds appear as rest points in 
McGehee coordinates just because the time transformation (2.2) is singular when E[h,r,a) = 0. 
Their counterparts in Cartesian coordinates correspond to turning points where qia) = 

5 Global dynamics in McGehee coordinates 

The typical orbit of equations (23) experiences a binary collision within a finite time, as we shall 
prove in this section. The exceptions constitute a set of zero Lebesgue measure. Moreover, there 
are no unbounded orbits, even when one introduces generalized solutions that allow for an orbit to 
pass through a binary collision. 
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5.1 Non-existence of unbounded non-colliding trajectories 

Let us start by ruling out unbounded, coUisionless orbits. In our problem orbits escaping to infinity 
cannot be ruled out with a simple argument based on the conservation of total energy h, because for 
any given h the distance r from the center of vorticity may grow without bounds while maintaining 
the same positive value of the (rescaled) kinetic energy i?, as should be apparent from Figures 3 
and 2. 

Definition 5.1 We say that a solution 7(C) := (''(C); ck(C)7 "0(0) of the system (23) is unbounded 

lim r{C,) ~ +00 

C-)- + oo 

Theorem 5.2 Equations (23) do not allow for unbounded coUisionless orbit. 

Proof. Let us consider the system (23); eliminating the time variable C this system is equivalent to 

r^ + 2 



(26) 



da 
dr 

d^ 
dr 



g tan('i/' — a) 

2 cos(V' + a) 



r 

r2 + 2 



rE 



sin(2a) cos{ip — a) 



tan^tp — a) 



By contradiction, we assume that there exists an unbounded motion r 1— > (a(r), ip{r)) starting from 
the initial condition ho = (rp, ap, "iAo) a-nd we suppose that the orthogonal projection of ho onto the 
(r, a)-plane is a point Hq lying in the region between the straight lines a ^ 0, r = r^ and the zero 
set of the function E, where r* is such that Eo{r^,,TT/A) = (see Figure 2). The case where ho lies 
between the straight lines a = 7r/2, r ~ r^ and the zero set of the function E is analogous and will 
not be explicitly addressed. In order to maintain i? > 0, an unbounded motion must satisfy 



(27) < a{r) < g{r), where g{r) 

for Eo{r) — 6{l — r^ logr). We observe that 



— arcsm 
2 



lim Eo{r)/2r^ 

r— >--[-oo 



lim log r 



4 cxp 



-00, 



Mr) 
2r2 



from which it follows that g{r) ^+00 l/(8r^) +o(l/r^). Thus r n- a{r) decreases to zero at infinity 
while keeping r 1— g(r) as an upper bound. From the first equation in (26) we have 



da 
dr 



+ 2 tan(V'(r) — a(r)) 



-\-oo 



and therefore 



lim 



tan(V-'(r) — a{r)) 



tan{'ip{r) — a{r)) 



0, 



We conclude that the function r 1— > a{r) should satisfy the limit equation 

da tan(?/'(r)) 
dr r 



From the second equation in (26) we have 



dr 



rE 



2 cos(?/' + a) 
sin(2a) cos('0 — a) 



tan(i/' — a) 



From the upper bound r t-^ a{r) G 0{l/r^) we also have 



r 






tan 


-^+00--^ 


a 


2 logr). 
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Therefore an unbounded motion should satisfy the limit problem 

tani/; 



3) 



The first equation is equivalent to 



da 
dr 

dip _ 
dr 6r log r 



1 



1 



ip{r) = — arctan(ra') 



hence unbounded orbits should asymptotically satisfy the following second order differential equa- 
tion 



(28) 



1 



1 + r^a'^ 6r log r 



1 + ra a' 



Multiplying the equation by a and taking into account that both a and a' go to zero as r — ^ +00 
it follows that the above equation reduces to 

1 



6r^ logr 



Observing that a > 0, it follows that an unbounded solution should be concave, which contradicts 
the inequality (27). □ 




Figure 7: Projection of the vector field in equation (23) on the (a, ip) plane at three different values 
of r. The black square and circle are, respectively the intersection with the plane of the ejection and 
collision heterocline of Proposition 4.4. The yellow vertical lines show the binary-collision planes 
at a = and a = 7r/2. The red and green lines are, respectively, the places where da/dC, = and 
dt/j/d^ = 0. The vertical grey lines show the intersection of the zero velocity manifold E = with 
the plane r = 1.25. The vector field is not shown in the dynamically meaningless region between 
the two lines, where E < Q. The regions shaded in light grey are those where dr/dC, < 0. Elsewhere 
it is dr/dC, > 0, except along the black dashed lines, where dr/dC = 0. 
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5.2 Generality of binary collisions 

Theorem 5.3 Any orbit not fully contained in the total collision manifold or in the zero velocity 
manifold, or in the central manifold of the rest point {r,a,ip) = (0, 7r/4, 57r/4) or in the stable 
manifold of the rest points S^-^ experiences a binary collision. Moreover, the union of the orbits that 
do not experience a binary collision is a set of zero Lebesgue measure in R'^. 

Proof. The total collision manifold and the zero velocity manifold are both two-dimensional, hence 
of zero Lebesgue measure. They are invariant manifolds and the orbits lying on them do not 
experience binary collisions, as illustrated in the previous section. 

Of all the points on the total collision manifold, only (r, «,■!/;) = (0, 7r/4, 57r/4) may be reached 
starting from r > (Remark 4.6). The heterocline 7724 belongs to the central manifold of (r, a, ip) = 
(0, 7r/4, 57r/4) associated to the eigenvector V2 of ^2 (Lemma 3.2). Owing to the non-uniqueness of 
the central manifold, other orbits are possible that become tangent to the heterocline 7724 as r — 0. 
Defining = '0 — a — tt, ^ = a — 7r/4, and approximating equations (23) at first order in /3, 0, the 
dynamics of orbits arbitrarily close to the heterocline is described by 

dr ^ r^ 



(29) 



dC r2 + 2 

dC ^ 

^ = -^•^(30 + 4/3). 

where we have omitted the irrelevant tt^^/^ factors. Here (/3, (j)) = (0, 0) is the heterocline, and r 
is a monotonically decreasing function of C. At each fixed value r > this system of equations 
has a contracting and an expanding direction on the {/3, (/))-plane, associated, respectively, to the 
eigenvalues Ai = — r(\/ 16E + Qr^ + 3r)/2 and A2 = r{\/l6E + 9r^ — 3r)/2. The contracting 

direction is given by the vector (r, /3, 0) = (0, 2E, r(\/ 16E + 9r'^+3r)), therefore the center manifold 
is a two-dimensional surface that becomes tangent to (r, a, ijj) = (r, a^a + n) at (a, ip) = (7r/4, 57r/4) 
as r — 7> 0. 

Another set of orbits that do not experience binary collisions are the stable manifolds of the 
points of ^3 associated to the stable eigenvector 1x2 (Lemma 3.3). These points are normally 
hyperbolic, and the union of their stable manifolds is, again, a smooth two dimensional surface. 
See [HiPuSh77] for details. 

Finally, from equations (23), in the rectangle R = (0, 7r/2) x [0, 27r] of the (a, ?/')-plane, we have 
the following properties of the orbits (see also Figure 7) : 

• a increases (with respect to time C) in the region of R where 

a < Tp < a + TT 

and decreases otherwise; 

• ijj increases (with respect to time C) in the region of R where 

arctan(TO(Q;)) < ip < arctan(m(Q;)) -|- tt. 

and decreases otherwise. 

• r decreases (with respect to time C) in the region of R where 

a + 7r/2 < V < a + 37r/2 

and decreases otherwise. 

The qualitative behavior of the vector field then rules out limit cycles or other invariant structures. 
Theorem 5.2 shows that collisionless, unbounded orbits are impossible. It then follows that orbits 
not on the above mentioned invariant manifolds, whose union is a set of zero Lebesgue measure, 
must reach a binary collision. □ 
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5.3 Generalized solutions: a boundedness result 



As a consequence of the fact that almost every orbit experiences a binary colhsion, we are lead to 
ask if it is possible to continue the solution after a binary collision in some natural manner. This 
is classical question in the context of n-body gravitational problems and it goes under the name of 
regularization of collisions. Motivated by the recent paper [CaTe] we shall introduce the notion of 
generalized solutions which, very roughly, are the continuation of the classical solutions after the 
binary collision in a suitable manner. In this section all "collisions" should be intended as "binary 
collisions", thus excluding the case of total collapse. Before proceeding further, we firstly recall 
some facts proven in [CaTe] in the case of weak central forces. 

One-center problem with logarithmic potential 

Let us consider the dynamical system associated with the conservative central weak force arising 
by a logarithmic potential having the singularity at the origin; i.e. the Cauchy problem: 



(30) 



q = dqU(q) 

(q(0),q(0))-(qo,Po)e(M2\{O}x 



where U{q) :~ — log(||<7||). A classical solution does not cross the singularity of the force, i.e. is 
a path q G "^^(T, \ {0}) where T denotes the maximal interval of existence. For the classical 
n-body problem, De Giorgi in [DeG96] proposed in to consider the smoothing of the potential as a 
regularization technique. An exhaustive analysis in the case of homogeneous potentials, including 
the classical Keplerian potential, was performed by Bellettini, Fusco e Gronchi in [BaFuGr03]. 
Following the idea of De Giorgi, the authors in [CaTe] removed the singularity at q = by smoothing 
the potential as follows: 

Uiq;e) = log ^\\q\\^+e^ 
and by considering the regularized problem 

q = dqU{q;e) 



(31) 



(g(0),q(0)) = (go,Po)e(R'xM2) 



Unlike in (30) the differential equation in (31) is no longer singular, so the initial value problem 
admits a global smooth solution q G "^""(R, M^) for every choice of the initial value (qo,Po)- Let 
Bq (R) be a ball of radius R centered at the origin of the configuration space, and let C x 
be the set of initial conditions within the ball leading to collision for the problem with e = 0. For 
every D G let qp G "^^(T, M^) be the collision solution where T denotes the maximal interval of 
existence. Denoting with q^.i/ the solution with initial data i/ and studying the asymptotic limit 
{s, v) — > (0, P), the authors in [CaTe] introduced the following two notions of regularization. 

Definition 5.4 We say that the problem (30) is weakly regularizable via smoothing of the potential 
in the ball Bq{R) if, for every v G -5^ there exist two sequences {sk)ki{vk)k tending to and v 
respectively such that there exists 



and the flow 



lim qe.^u, = q 



q{t) vG,9' 



Quit) . 

is continuous with respect to v (i.e. the initial point) 



Definition 5.5 The singular one center problem (31) is strongly regularizable via smoothing of the 
potential if there exists R such that for every i? G the limit 

lim q^^i, = q 
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q(t) v(^.9' 



exists and the flow 

is continuous with respect to v. 

The authors in [CaTe] then prove that 

Proposition 5.6 The logarithmic one center problem is globally regularizable via smoothing the 
potentials according to Definition 5.5. 

As a consequence of this result, it follows that it is possible to continue the solutions after the 
collision by transmission, in the following sense. 

Definition 5.7 Let t i— > qo(t,), t € [0,ro) be a collision path, and Tq the collision instant. Define 
the transmission solution go for t G [0, 2Tq\ as follows: 

q{t)=qo{t) <e[0,To] 

q{t) ^ ~qo{2TQ - t) te[To,2To] 

This is the right definition to set in order to have the continuity with respect to the initial data of 
the flow obtained by replacing the collision solution go with the transmission solution. 

Local regularization of collisions and generalized solutions 

We apply the above results to the solutions of our problem that experience binary collisions. The 
presence of other vortices, of course, complicates the picture, but it should be intuitively clear that 
any binary collision, locally, is just a central problem, as we show formally below. 

In Cartesian coordinates, let us consider instead of the singular potential defined in (13), the 
non-singular one given by 



(32) U{qo-e) = - log (Vll^o -.99011'+. 

that substitutes (13) in the equations of motion (5). For the Klein group D2 (see section 3), the 
above expression reduces to 

(33) C/(<7i, 92; e) = \ log ((4q? + e^){Aql + e^){A \\q4^ + e^)) . 

We recall that a simultaneously binary collision in our problem corresponds to a solution q of the 
Newton's equations (4) 

rq = — 

dq 

for U{q) := — log(8||<7||(7i (72) such that at some time instant cr* one and only one coordinate of the 
point q{(J*) (91 (o"*), (72(0'*)) is zero. This means that the support of the solution intersects one 
of the coordinate axis qi or (j2. Without loss of generality we can assume that at some instant 
we have (72(0'*) = 0. Thus there exists 6 := 6{(T^,) > such that 

- log((j2 A'l) < U{q) < - log(<72 K2), 

where Ki := 8 and K2 := 8 (C— 5)^ and C := gi(cr*). Therefore, locally in the neighborhood 

of a singularity, the potential is, up to a constant, the central logarithmic potential described in the 
paragraph above. 

Remark 5.8 By continuing the solution after a binary collision as described above we can replace 
the colliding trajectory, which does not exists beyond tr*, with a transmission solution that exists 
also at later times. However at some instant this new trajectory could experience another binary 
collision that can, again, be extended by transmission, and so on for any number of collisions. 
Although we could construct arbitrarily long solutions containing an infinite number of collisions, 
the above arguments do not guarantee the uniform convergence of the generalized solutions. The 
extension of a generalized solution up to an infinite time ( if possible ) would require a much 
and careful analysis than that presented above. 
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Definition 5.9 A generalized solution q € "^^(J) n '^'^{J) of the problem is a classical solution 
("classical" meaning a curve that satisfies the equations (4) pointwise) on the time interval 
J = J\Y , where Y :~ {cti, . . . .a„} is a finite set of times. At those times the solution q experiences 
a binary collision and is continued by transmission as explained above. 

Before proceeding further we translate the notion of generalized solution to the McGehee setting. 
In these coordinates, a solution experiences a binary collision when it crosses either the plane a = 
or the plane a = 7r/2, corresponding, respectively, to collisions with 52(0'*) = and gi(cr*) = in 
Cartesian coordinates. The graph of r 1-^ a{r) on the (r, a) plane for a transmission orbit colliding 
at r = is locally symmetric with respect to r*. The graph of r i— > ip{r) has a jump discontinuity 
of TT at r = r«. 



Boundedness of the generalized orbits 

The equations of motion given in (23) may be written as follows: 

r-2 + 2 



(34) 



da 
dr 

d4 
dr 

<hp_ 
da 



tan(?/' — a) 



j1±1 

rE 

E 



2cos(V' + a) 
sin(2a) cos('0 — a) 

2 cos{'ip + a) 



tani^l) — a) 



sin(2a) sin('!/) 



- 1 



As observed in Theorem 5.2, if r 



-00 then necessarily a — )■ 0, otherwise the constraint E > 



would be violated. In the following we shall distinguish between two cases: 

1. if a{r) G 0(l/r^) for r — >■ 00 then Eo{r,a) w c for some positive constant c; 

2. if Q!(r) e o{l/r^) for r ^ 00 then EQ{r,a) w — 2r^loga. 
First case. For fixed we have that 

tan?/! 



tSLn^ip — a) 

cos{'ip — a) = cos tp + a sin ip + o{a), 
and hence the system given in (34) reduces to 



■ + o(a), 

1 + a tan ip 

cos(V' + a) = costp — asin^p + o(a). 



(35) 



da 
dr 

d4_ 
dr 

dip 
da 



- tan(-!/') 
r 



2 tanip 
- 2a tan ip 



atanip 



Lemma 5.10 Let us consider an orbit starting at the point Pq = (ro,ao, V'o) o.'f^d let us also assume 
that tan?/; is bounded. Then the projection of the orbit onto the {a,ip)-plane intersects the line 
a = for r > ro if tan -0 < and for r < r^ otherwise. 

Proof. The result immediately follows from the first equation in (35). □ 
In order to show that generalized solutions are bounded, we must consider the following three cases: 



(i) If atanip e 0(1) then the system (35) reduces to 
(36) 



da 
dr 



1 

a r ' 



dr 



dp 
da 
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(ii) If a tamp g o(l) then wc have 



(37) 



da 
dr 



taTiip 



dr 



dtp 

da a tan 



(iii) If a tan -tp is unbounded then 



(38) 



da 
dr 



tanip 



dip_ 
dr 



rtanip, 



(hp_ 
da 



r\ 



An interesting consequence of Lemma 5.10 is that a gencrahzed solution may temporarily exhibit 
a monotonicity (with respect to r) of the collision points. More precisely, taking into account the 
above estimates for dip/dr, the following result holds. 

Lemma 5.11 Under the assumptions of lemma (5.10) it follows that the projection of the orbit onto 
the {a,^)-plane intersects the line a = at the instants (r,;)igiv with ri < r,;+i if tan ^ < and 
ri > r^+i otherwise. 

However, the same asymptotic estimates above imply that a monotonic sequence of collisions cannot 
be arbitrarily long. Otherwise stated, tan^j is not bounded along the whole orbit. 

Lemma 5.12 If a(r) G 0(1/?''^) unbounded generalized solutions do not exist. 

Proof. By contradiction, we observe that in an unbounded solution 7 the first component r i— > a(r) 
tends either to or to 7r/2 for r — )- 00. Without loss of generality we discuss only the case a = 0. 
We also observe that along an unbounded generalized solution ^ — > when r — > +00. We are 
also assuming that a{r) £ 0{l/r^). This immediately rules out case (i) above, as it would imply 
^ +00. In case (ii) ip is unbounded because of a{r) G 0{l/r^). For the same reason, in case 
(iii) ta'n{ip) is unbounded. Therefore there exists a r, such that ipi'''*) = '''/2 + kn for some fc € Z. 
For r > r, the function r i— > a{r) decreases and the monotonically growing sequence of collisions 
stops. □ 
Second case. In this case the system (34) reduces to 

da 1 , 
— ^ - tan(V') 
dr r 



(39) 



di; _ 1 
dr 2r log a 

dip _ 1 
, da log a 



2tan'i/; 



1 — 2a tan -0 
atan-i/) 



By arguing exactly as in the previous case, in this case is also possible to show that the following 
result holds. 



Lemma 5.13 If a(r) G o(l/r'^) unbounded generalized solutions do not exist. 



6 Global dynamics in physical coordinates 

From now own, when we refer to solutions we mean generalized solutions that may go through 
binary collisions. 

The transformations linking the physical coordinates (gi, (72) and momenta {pi,P2) to the McGe- 
hce coordinates (r, a, ip), for a given value of the total energy h, are 



(40) 



(41) 



qi = r e ' cos a 
q2 = r e ' sma 



pi = zi/r ~ y Eir/r cos ip 
P2 ~ Z2/r = \/ Eir/r sinip 
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We recall that the various time scalings introduced along calculations have the effect that both the 
invariant manifold corresponding to the total collapse and that corresponding to zero velocity are 
reached asymptotically as the rescaled time goes to infinity. 

We further observe that not all the rest points on the total collision manifold have a physical 
meaning. Only the points {r,a,ip) = (0, 7r/4, 7r/4) and {r,a,ip) = (0, 7r/4, 57r/4) are physically 
relevant. They correspond to central configurations in which each vortex lies on the vertex of a 
square, having the barycenter at the origin and the vertices on the bisectrices of the coordinate 
axis. Thus the two heteroclinic connections of Proposition 4.4 joining the total collision manifold 
and the zero velocity manifold correspond, respectively to the homographic ejection orbit from total 
collapse to the zero velocity manifolds and to the homographic collision from the zero velocity to 
the total collapse manifolds. 

This is illustrated in Figure 8 by a numerical solution of the equations of motions in physical 
coordinates using the non-singular potential (33) for an initial condition close to the ejection hete- 
rocline. Initially the solution follows closely the heteroclinic cycle shown in figure 6. Eventually, the 
solution leaves the collision heterocline before reaching the total collapse (which is a saddle point) 
and it is subject to a sequence of two binary collisions close to total collapse. 

In Figure 9 we illustrate a case complementary to the above one. We show a numerical ap- 
proximation to an orbit lying in the center manifold of the total collision point at {r,a,Tp) = 
(0, 7r/4, 57r/4). The orbit approaches the homothetic configuration (which, in McGehee coordinates 
is the heteroclinic connection 7724) with a projection on the (a, ip) plane which is tangential to the 
line of rest points J3^2- In Cartesian coordinates our choice of initial conditions corresponds to 
the four vortices initially close in pairs, as if they were just emerging from a binary collision, but 
with the special choice of the momenta that sends them on a trajectory tangent to the homothetic 
solution. 

Summarizing we have 

Theorem 6.1 The only homothetic ejection/ collision orbit from/to total collapse is that homothetic 
to the planar central configuration, which is bounded. 

Proof. The existence and boundedness of the ejection/collision orbit homothetic to the planar 
central configuration is a direct consequence of Proposition 4.4. For a hypothetical solution where 
the position of the vortices is at all times a homothetic transformation of the vertices of a given 
rectangle, then the momentum of each vortex would be parallel to the line joining the vortex and 
the origin. However, it is clear that only at a = 7r/4 all the forces are balanced, therefore we exclude 
other solutions. □ 

Remark 6.2 Recalling from Sundman type estimates proven in [BaFeTeOS] that, of all the points 
on the collision manifold, those corresponding to a homothetic ejection/collision (namely {a,ijj) = 
(7r/4, 7r/4) and {a^ijj) = (7r/4, 57r/4) j are dynamically reachable from outside the manifold, then 
the other rest points on the collision manifold and the heteroclinic connections between them are 
dynamically meaningless. 

Theorem 6.3 Each solution different from the homothetic central configuration experiences a bi- 
nary collision in finite time. Moreover, any solution starting arbitrarily close to the homothetic 
ejection solution has an orbit that may not reach the upper bound of the homothetic orbit before 
experiencing a binary collision. 

Proof. This is a direct consequence of the theorem 5.3. □ 
When the four vortices are arranged at the vertex of a rectangle with a large ratio of side 
lengths, then there are solutions in which the pairs of vortices undergo repeated binary collisions 
along the same axis. In particular, for some of these solutions the sequence of collisions initially 
moves away from the origin of the Cartesian axes. However, unbounded solutions are not possible, 
and eventually the sequence reaches a turning point and returns towards the origin. An example 
of this dynamical behavior is shown in figure 10. When seen in the McGehee coordinates, these 
are solutions advancing (up to a turning point) into the narrow funnel between the zero velocity 
manifolds of adjacent domains. In figure 11 we show a numerical solution where a sequence of 
binary collisions approaching the origin from the right undergoes a single binary collision along the 
q2 axis and the continues moving leftward on the negative qi semi-axis, although with a markedly 
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Figure 8: A near-ejection orbit starting at (r, a, V') — (0.3, 7r/4 + 10~^, 7r/4) with h ~ 0. The initial 
conditions are transformed to physical coordinates using equations (40) and (41), then a numerical 
integration of the equations (5) with the potential (33) and e = 10~^ generates the orbit. Panel A): 
the orbit in physical coordinates, projected on the (91,(72) plane. The inset magnifies the region 
close to the origin (corresponding to the total collapse). Panels B), C), D): the orbit transformed 
back to McGehee coordinates, projected respectively on the {r,a), {a,ilj), and {r,ip) planes. For 
clarity, after each binary collision the orbit is extended with continuity into the adjacent domain, 
rather than into the fundamental one, which would introduce visually bothersome jumps in the ip 
coordinate. The circle shows the position of the initial condition. The yellow lines at a = 0, tt, • • • 
correspond to binary collisions at qi ~ 0, those at a ~ 7r/2, 37r/2, • • • to binary collisions at q2 ~ 0. 
The grey lines in panel B) are the points with E = 0. 
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Figure 9: An orbit as close as numerieally possible to one lying on the center manifold of the 
total collision point (r^a^ip) = (0, 7r/4, 57r/4). The orbit starts approximately at (r, a,";/") = 
(1.00264534 • • • , 0.04149204 • • • , 2.54344846 • • •]) with /i = 0. The dashed lines are the homoth- 
etic orbit corresponding to the heterocline 7724- The red and the green curve in panel C) are, 
respectively, the projection on the (a, V') plane of the lines of rest points 3^2 and Details of 
the numerical solution and panel description as in figure 8. 

different amplitude. Eventually, this it will reach a turning point, and return towards the origin. 
Solutions that oscillate along an axis bouncing back and forth between turning points of opposite 
sign are not uncommon. However, we were unable to maintain the oscillations, which are aperiodic, 
for arbitrarily long times. Eventually the orbit spends some time winding closely around the origin 
at distances corresponding to values of r such that r, a) > for any a. After these complicated 
and chaotic-looking transients the orbit may resume its oscillations along one (not necessarily the 
same) axis. 

We may formalize these findings as follows. 

Theorem 6.4 Ij the initial configuration is a rectangle with large enough aspect ratio, the configu- 
ration evolves by passing through an arbitrary number of simultaneous binary collisions between the 
pairs of closest vortices. 

Theorem 6.5 The projection of any solution in the configuration space is hounded. 

Proof. This is a direct consequence of theorem 5.2. □ 

Theorem 6.6 The set of initial conditions leading to total collapse or leading in an infinite amount 
of time to the outermost boundary region where the motion is allowed has zero Lebesgue measure. 

From the physical point of view the outermost boundary region among the others depends on the 
total energy h of the system and it is given by the hypersurface having the property that any 
solution of the Newton equations after reaching with zero velocity, the solution orbit falls down. 

Finally, we may go back to the initial interpretation of the equations of motion (4) as steady so- 
lutions of the equations (1) for infinitely tall, nearly parallel, vortex filaments in a three-dimensional 
space. In that setting our time-like variable really is the arc-length parameter along the vortex fil- 
ament. From this point of view an orbit of our dynamical system describes the mutual position of 
four steady vortex filaments at various heights. The points of intersection of the vortex filaments 



25 





26 



with an horizontal plane, and the inclination of tangent to the filaments at the intersection points 
are the parameters that completely determine the shape of the whole filament. In particular, some 
of our findings about the equivariant solutions subject to the D2 symmetry may be recast as follows: 

1. If the vortex filaments join together into a total collapse, or separate from collapse, then they 
do so by approaching a square configuration; 

2. the distance among the vortex filaments is bounded at all heights; 

3. binary collisions of the vortex are generic, that is, the set of intersection points and tangents 
generating non-colliding filaments has zero Lebesgue measure. 

7 Further perspectives and closing remarks 

Among all the solutions of a system of vortex filaments some special and very important class is 
represented by the so-called helicoidal solutions, namely solutions of (1) of the form. 

for some complex valued function q : M ^ C. In the special case in which q{a) ~ A(o')^ it is called 
a homographic solution and if if X{a) = e"^'' we shall refer as a relative equilibrium. We observe 
that, for fixed frequencies w, v the function 

only depends on the choice of the configuration ^. Moreover every filament has the same shape of 
helicoidal type and the time evolution corresponds to a rotation equal for each vortex filament of 
the initial configuration. Therefore the solution ^' is periodic in a and constant in time. 

1. An interesting question that should be addressed is to consider helicoidal solutions instead 
of stationary solution. In this case the dihedral potential we need to add a quadratic term. 
However we think that our techniques should work also in this case. 

2. Of some interest in the applications is the study of a singular logarithmic type potential on a 
more general surface without boundary. 

As already observed in the introduction this transformation was firstly introduced in Celestial 
Mechanics by R. McGehee in its celebrated paper [McG74]. Furthermore this technique became a 
milestone in order to investigate the orbital structure of the n-body problem in the neighborhood 
of total collapse. In fact it was employed in several problem giving a lot of important feature on 
the global dynamics of this interesting problem. However in that context, due to the fact that the 
potential is a homogeneous function (of degree —1) allow us to decouple the system in a scalar 
equation containing the radial part and in a system which take into account the angular part. 
Therefore it makes this technique more feasible for the application. In fact it is possible to study 
separately the angular and radial part; in this perspective it is possible to read the full behavior 
of the system by looking only at the angular part. All of this breakdown in our context since the 
potential is not homogeneous anymore. 

Another important difference which makes things more involved is that the equilibrium points 
are not isolated (in fact they appear in family) and they do not have a hyperbolic character (and 
only some of them are normally hyperbolic manifolds). For instance, to the knowledge of the 
authors it is not known if there exists a sort of Palis Inclination Lemma in this context. 

However since the differential of the potential is a homogeneous function, this still guaranteed 
the existence of self-similar (homothetical, homographical) motions. 

As direct consequence of the energy relation, for each fixed energy level h the motion is confined 
between the two invariant manifold (£" = 0, r = 0) and the hyperplanes corresponding to the 
simultaneously binary collision. However as already observed the region in which the motion is 
allowed is not compact. This is completely different for instance with the case studied by the 
authors in [StFoOS] and in our case a priori unbounded motion can happen. However as proved 
above there are no unbounded orbits. 

An intriguing question is to understand if this weak logarithmic singularities always prevents 
the existence of unbounded motions. 
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